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Multivariate C0GARCH(1, 1) processes are introduced as a continuous-time models for mul- 
tidimensional heteroskedastic observations. Our model is driven by a single multivariate Levy 
process and the latent time- varying covariance matrix is directly specified as a stochastic process 
in the positive semidefinite matrices. 

After defining the COGARCH(l, 1) process, we analyze its probabilistic properties. We show 
a sufficient condition for the existence of a stationary distribution for the stochastic covariance 
matrix process and present criteria ensuring the finiteness of moments. Under certain natural 
assumptions on the moments of the driving Levy process, explicit expressions for the first and 
second-order moments and (asymptotic) second-order stationarity of the covariance matrix pro- 
cess are obtained. Furthermore, we study the stationarity and second-order structure of the 
increments of the multivariate COGARCH(l, 1) process and their "squares". 
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1. Introduction 

In this paper, a multivariate extension of the continuous-time generalized autoregressive 
conditional hcteroskedasticity (COGARCH, for short) process of order (1, 1) introduced 
in [26] is defined and studied in detail. The one-dimensional C0GARCH(1,1) process 
(see also [10, 22, 27]) is given as the solution of 

AGt = ^dLt, (1.1) 
dvt = -P{vt- -c)dt + avt- d[L,i]?, (1.2) 

using the discontinuous part [i, L]'' of the quadratic variation of a univariate Levy process 
L, parameters a, /?, c > and initial values Go = 0, > 0. The process G is referred to as 
the C0GARCH(1,1) process and the variance process v as its volatility process^ where 
the name "volatility process" derives from the term typically used in economics. 
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Heteroskedastic data are often modelled with (normal) variance mixture models. In 
such a model, one has X„ — ^/v^En for n G N, where e is an i.i.d. sequence and v a 
sequence of positive random variables modelling the current variance of the observations 
X. Typically, one has, moreover, that £„ and u„ are, for fixed n S N, independent. Ob- 
viously, equations (1.1) and (1.2) constitute a continuous-time counterpart of a variance 
mixture model with some special process for the variance and driven by a single Levy 
process. Loosely speaking, the increments dLt are "mixed" with the variance vt- and 
the two are independent for fixed t. 

In a multivariate setup, the positive variance v needs to be replaced by a covariance 
matrix process V. Thus, the volatility process has to be a stochastic process in the pos- 
itive semidefinite matrices. This requirement leads to challenging questions in modelling 
and interesting mathematical issues since, in particular, only very few continuous-time 
stochastic processes in the positive semidefinite matrices have been thus far studied 
(mainly various "Wishart" processes, see [11, 12, 14, 19, 20]; or, recently, Ornstein- 
Uhlenbeck-type processes, see [5, 34, 35]). Appropriate multivariate models for het- 
eroskedastic data are, however, clearly needed because in many areas of application, 
one has to model and understand the joint behavior of several time series exhibiting non- 
trivial interdepcndencies. Moreover, for various reasons (for example, unequally spaced 
observations, inference at several frequencies or amenability to continuous-time financial 
theory), it is often desirable to use continuous-time models instead of related discrete- 
time models like GARCH models, for instance. 

After briefiy stating some preliminaries regarding notation and Levy processes in Sec- 
tion 2, we introduce our multivariate C0GARCH(1,1) processes (MUC0GARCII(1, 1), 
for short) in Section 3 and establish well-definedness. Thereafter, we analyze its volatil- 
ity process in Section 4. In the first part of that section, we present a univariate 
C0GARCH(1,1) process that bounds the volatility process in a norm intrinsically re- 
lated to the autoregressive parameter and use this bound to give sufficient conditions for 
the finitcness of moments. This is followed by a demonstration that the volatility process 
alone and the MUC0GARCH(1, 1) process together with its volatility are strong Markov 
processes. Moreover, we establish conditions for the existence of a stationary distribution 
of the volatility in Section 4.2. In the last part of Section 4, we calculate the second-order 
structure of the volatility process explicitly under certain assumptions on the moments 
of the driving Levy process and establish (asymptotic) second-order stationarity. 

In Section 5, we focus on the increments of the MUG OG ARCH ( 1, 1) process itself, 
showing that it has stationary increments provided the volatility is stationary. Thereafter, 
we calculate the second-order moment structure of the increments (that is, the returns 
in a financial context) observed on a regularly spaced discrete grid and their "squares" 
(that is, the increments times their transposes). Here, we obtain, in particular, that the 
increments have zero autocorrelation, but their "squares" have exponentially decaying 
autocorrelation. Moreover, the explicit expressions for the moments obtained make the 
processes amenable to statistical estimation. 

Finally, we present all proofs, together with auxiliary technical results, in Section 6. 
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2. Preliminaries 

2.1. Notation 

We denote the set of real mx n matrices by Mm,n{^)- If m = n, we simply write M„(M) 
and denote the group of invcrtiblc n x n matrices by GL„(M), the linear subspace of 
symmetric matrices by §„, the (closed) positive semidefinite cone by and the open 
positive definite cone by §^^. /„ stands for the n x n identity matrix. The natural 
ordering on the symmetric n x n matrices shall be denoted by <, that is, for A,B £ 
Af„(R), we have that A<B\i and only if B - ^ e S+ (likewise, A<B^B-A& §++). 
The tensor (Kronccker) product of two matrices A,B is written as A®B. vcc denotes the 

2 

well-known vectorization operator that maps the set of n x n matrices to R" by stacking 
the columns of the matrices below one another. For more information regarding the tensor 
product and vec operator, we refer to [24], Chapter 4. The spectrum of a matrix is denoted 
by a{-) and the spectral radius by p(-). Finally, A* denotes the transpose (adjoint) of a 
matrix A £ Mm,n(^)- 

Norms of vectors or matrices are denoted by || • || . If the norm is not further specified, 
then it is irrelevant which particular norm is used. 

Throughout, we assume that all random variables and processes are defined on a given 
filtered probability space (fi, J", P, (J^t)tGr)7 with T = N in the discrete-time case and 
T = in the continuous-time one. Moreover, in the continuous-time setting, we assume 
the usual conditions (complete, right-continuous filtration) to be satisfied. 

Furthermore, we employ an intuitive notation with respect to (stochastic) integra- 
tion with matrix-valued integrators, referring to any of the standard texts (for exam- 
ple, [36]) for a comprehensive treatment of the theory of stochastic integration. Let 
(^t)teR+ Mjn,n{^) and {Bt)teR+ -^r,s(R) be cadlag and adapted processes and 
{Lt)teR+ ill Mn^rO^) be a semimartingale. We then denote by J^^ Ag- dLgBg- the ma- 
trix Ct in Mjn,sO^) which has ijth element Cy,t = Y.k=iY^'i=i lo ^ik.s-Bij^s- dLki,s- 
Equivalently, such an integral can be understood in the sense of [32, 33] by identify- 
ing it with the integral As- dLg, with At being, for each fixed t, the linear operator 
M"„_^(R) ^- M™,s(M),X H> AtXBt. If {Xt)teR+ is a semimartingale in R"* and {Yt)t(zR+ 
one in R", then the quadratic variation {[X,Y]t)teM+ is defined as the finite variation 
process in Mm,n(R) with components [X,Y]ij^t = [^iiYjjt for t € R+ and i = 1, . . . ,m, 
j = l,...,n. 

2.2. Levy processes 

Later, we shall use Levy processes (see [1, 36, 39], for instance) both in M.'^ and in the 
symmetric matrices Sd- Thus, we briefly recall the relevant basic notions for them now. 

We consider a Levy process L = {Lt)t(rR+ (where Lq = a.s.) in R'' which is determined 
by its characteristic function in the Levy-Khintchine form iJ[e'^"'^*^] = exp{i-0i(u)} for 
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where 7^ G M'', tl € and the Levy measure vl is a measure on R'' satisfying 
z^l({0}) = and /jgd(||a;p A 1)i'l{'^x) < 00. Moreover, (•,•) denotes the usual Euchdean 
scalar product on Mf^. 

We always assume L to be cadlag and denote its jump measure by /xl, that is, fiL is the 
Poisson random measure on M+ x M'' \ {0} given by hl{B) = tt{s > : (s, — Lg-) S B} 
for any measurable set B C K+ x R'' \ {0}. Likewise, /ii(ds, da;) ~ /XL(ds, dx) — ds i'L{dx) 
denotes the compensated jump measure. 

If ^ then L has finite mean and covariance matrix given by 



Provided = and /||3,||<i ll-'^ll'^L(da;) < 00, the Levy process L has paths of finite 
variation and 



Regarding matrix-valued Levy processes, we will only encounter matrix subordinators 
(see [3]), that is. Levy processes with paths in Sj. Since matrix subordinators are of finite 
variation and tr{X*Y) (with X,Y GSd and tr denoting the usual trace functional) defines 
a scalar product on Sd linked to the Euclidean scalar product on R'^ via ti{X*Y) = 
vec(X)* vec(F) = (vec(y), vec(X)), the characteristic function of a matrix subordinator 
can be represented as 



with drift 7^ S §J and Levy measure u^. 

The discontinuous part of the quadratic variation of any Levy process L in R*^, 




(2.1) 




£;(eit'-(i*^))=exp(tVL(^)) 





is a matrix subordinator with drift zero and Levy measure given by 




for all Borel sets S C S^. 
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3. Definition of multivariate COGARCH(l,l) 
processes 

The main idea for the definition of a multivariate C0GARCH(1, 1) process is to replace 
the noise e of a multivariate GARCH(1,1) process (see [6, 16] and references therein) 
by the jumps of a multivariate Levy process L, and the autoregressive structure of the 
covariance matrix process by a continuous-time autoregressive (AR) structure (that is, 
an Ornstein-Uhlenbeck (OU) type structure). So, the idea is again basically the same as 
in [10] for the univariate COGARCH(p, q) process. 

In the simplest BEKK GARCH(1, 1) model of [16], the volatility process is given by 

S„ = C + AI]y„V„-i4_iSy_'iA* + BI]„_iB* (3.1) 

with C G Sj^, A,B^ and (en)neNo being an i.i.d. sequence in W^. This shows 

that the dynamics of (E„)„gii}o those of a multivariate AR process, which is 
"self-exciting" in the sense that we have an AR structure with the noise given by 

(Sy_;^e„_l£*_lE/_;^)„gN- 

Replacing the AR structure with an OU-type structure, using V^^J.^ d[L, L]^ V^/^ as 
"noise" , where L is a d-dimensional Levy process, and using the same linear operators 
as for positive scmidefinite processes of OU type (see [5, 35]) now leads to a multivariate 
continuous-time GARCH(1,1) process G (referred to as a MUC0GARCH(1, 1) process 
in the following) given by the following definition. 

Definition 3.1 ^UCOGARCH(l, 1);. Let L he an W^-valued Levy process ana 



A, Be Md{M.), C G §J+. The process G = {Gt)tm+ solving 

dGt = VlI^dLt, (3.2) 

Vt = G + Yt, (3.3) 

dYt = {BYt^ + Yt-B*) dt + AVll^ d[L, L]] V^l'^A*, (3.4) 



with initial values Gq in M'' and Yq in Sj", is then called a MUC0GARCH(1, 1) process. 

The process Y = (lt)tgR+ with paths in is referred to as a MUC0GARCH(1, 1) 
volatility process. 

As we are only dealing with MUCOGARCH processes of order (1, 1), we often write 
only "MUCOGARCH" instead of "MUC0GARCH(1, 1)" in the sequel. 

We can also directly state a stochastic differential equation (SDE) for the covariance 
matrix process V: 

dVt = {B{Vt^ -C) + {Vt- - C)B*) dt + AVI/'^ d[L, L]l Vtl'^A*. (3.5) 

This SDE has a "mean-reverting structure" (provided cr{B) C (—00,0) -l-iM), namely, V 
returns to the level G at an exponential rate determined by -B, as long as there are no 
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jumps. However, since all jumps are positive semidefinite, as we shall see, C is not a 
"mean" level, but a lower bound. 
Equivalently, wc can use the following representation using the vec operator: 

d vec(yf ) :=^{B®I + I®B) Ycc{Yt-)dt + {A® A){Vtl^ ® d vcc([i, i]^^), 

d vec(Vt) ^{B®I + I® B){vcc{Vt^ ) - vec(C)) dt + {A® A){vl!^ ® d vcc([i, L]^). 

For the MUCOGARCH process to be well defined, it is necessary that there exists a 
unique solution to the above system of stochastic differential equations and that V does 
not leave the set §J. In the following, it is implicitly understood that our processes and 
stochastic differential equations are not living on the space Md(M) (resp., M''^), but on 
the linear subspace Sd of symmetric matrices (resp., vec(Sd))- The latter can, as usual, be 
identified with M''(<^+i)/2^ when appropriate. The importance of this lies in the fact that 
Sj"*" is an open subset of §d and it is most natural to consider (stochastic) differential 
equations on open sets. 

Theorem 3.2. Let A,B £ Md{R), C e S^^ and L be a d- dimensional Levy process. 
The SDE (3.4) with initial value Yq € Sj' then has a unique positive semidefinite solution 
(^t)tGR+. The solution (Yt)t£M+ locally bounded and of finite variation. Moreover, it 
satisfies Yt > e^'Foe^'* for all tGM.+ . 

Remark 3. 3. (i) An analogous result holds when considering C G and restricting the 
initial value to Yo <E ( "locally bounded" needs to be replaced with "locally bounded 
within §J", as defined in [5], Definition 3.1). All of the following results (except those 
regarding the existence of stationary solutions) can be immediately adapted to this case. 

(ii) For d = 1, it is straightforward to see that our definition agrees (after a 
reparametrization) with the case p = q = 1 of the general COGARCH(p, q) definition 
given in [10], that is, Yt agrees with their process aiYt and our Vt with their Vt+. Hence, 
[10], Theorem 2.2, implies that our definition agrees with the original definition given in 
[26]. 

Likewise, we could have considered the SDE (3.5). Using the relationship between (3.5) 
and (3.4), we obtain the following. 

Corollary 3.4. Let A,B AId{M.), C G and L be a d-dimensional Levy process. 

Assume that the initial value satisfies Vq > C The SDE (3.5) then has a unique positive 
definite solution (Vt)teR+ and Vt>C + e^*(Vb - C)e^** > C for all < G K+. 

It may appear natural also to allow initial values Vq G '§(^ with Vq < C. In this case, 
one still has that Vt > C + e^*(Vb — C)e^ * needs to be true for any solution of (3.5), as 
long as it exists. However, in this case, the solution of (3.5) may leave the set Sj" and 
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thus have only a finite lifetime, as the following example shows. Take 

[2 0\ /0.5 

^-[02)^ ^°-[o 0.5 

/-0.51n(10/9) \ fl 

\ 1 -0.5^10/9)^ ' i^l 

and L as the zero Levy process. We then obtain that 

So, 1/1 ^ although Vb G Note that this problem also arises with positive probability 
if the driving Levy process is compound Poisson, as it may then happen that there is no 
jump until time 1. 

Remark 3. 5. The insight gained from positive semidefinite OU processes in [5] suggests 
that all eigenvalues of B should have negative real part if one wants a stationary COGA- 
RCH volatility process as covariance matrix process. It is clear that, in this case, Vt — ?> C 
as i — )■ 00 if the Levy process had no jumps. Thus, in general, the process V tends to 
C, as long as the driving Levy process does not jump. The above counterexample shows 
that when the process V is smaller than C (in the ordering of the positive semidefinite 
matrices), this does not occur in a "straight" manner, whereas in the univariate model, 
the volatility process is always increasing below C (see [27], Proposition 2). 

Similar to the usual shot noise representation of OU-type processes, we have the fol- 
lowing. 

Theorem 3.6. The MUC0GARCH(1, 1) volatility process Y satisfies 

Yt = e^'Foe^** + /* e^(*-^)A(C + Ys^)'/^ d[L,L]l {C + n_)i/2A*e^*(*-^) 

, ^° (3.6) 
for allteR+. 

Recently, [38] studied univariate equations of the form X{t) = J{t) + Jq g{t — 
s)f{Xs~)dZs and their relation to certain SDEs. In particular, they obtained uniqueness 
of the solutions under uniform Lipschitz assumptions on /. Our equation (3.6) is a mul- 
tivariate equation of this type, with / being only locally Lipschitz. From the arguments 
given in [38], one sees that their Theorem 5.2 remains valid in a multivariate setting. 
Using a localization procedure as in the proof of [40], Theorem 6.6.3, this uniqueness 
result extends to / being defined only on an open subset and locally Lipschitz. Hence, 
(3.6) provides an alternative characterization for the MUCOGARCH volatility process. 

So far, we have excluded the MUCOGARCH process G itself from the analysis. How- 
ever, the following result is obtained along the same lines as Theorem 3.2. 
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Theorem 3.7. Let A,B G Afd(R), C € S^"*" and L be a d- dimensional Levy process. The 
system of SDEs (3.2), (3.4) then has a unique solution {Gt,Yt)i.fz]s^+ with paths in M."^ x Sj" 
for any initial value {Go,Yo) in M."^ xS'^ . The solution (Gt, ^)tGR+ semimartingale 
and locally bounded. 



4. Properties of the volatility process 
4.1. Univariate COGARCH(l,l) bounds 

Wc now show that, similarly to the COGARCH(p, (7) case (see [10], Lemma 9.1), the 
norm of a MUC0GARCH(1, 1) volatility process can be bounded by a univariate 
C0GARCH(1,1) volatility process. This immediately gives useful conditions for the 
finiteness of moments and has far-reaching implications regarding the existence of sta- 
tionary distributions. 

In the following, we shall consider a special norm that fits our model particularly 
well. II • II2 denotes the operator norm on Af(j2(K) associated with the usual Euclidean 
norm. Assume, now, that B is diagonalizable and let S G GLd{C) be such that S~^BS is 
diagonal. We then define the norm || • \\b,s on Md2{R) by \\X\\b,s ■= WiS'^ S-^)X{S 
5") II 2 for X G Mj_2{M.). It should be noted that || • \\b,s depends both on B and on the 
choice of the matrix 5* diagonalizing B. Actually, || • ||s,s is again an operator norm, 
namely the one associated with the norm ||x||b,s := II ('S'"^ iS'~^)a;||2 on R'' . Besides, 
II • ||b,S is simply the norm |j • II2, provided S' is a unitary matrix (see [23], page 308). 

Theorem 4.1. Let Y be a MUCOGARCH volatility process with initial value Yq e §J 
and driven by a Levy process L in R''. Assume, further, that B £ M(j(M) is diagonalizable 
and let S G GLd{C) be such that S^^BS is diagonal. The process solving the SDE, 



dyt = 2\yt- At+\\S\\l\\S-^\\lK2,B\\A®A\\B^s[ ^ + Vt- ) dL*, 



yo = II vcc(yo) 



(4.1) 



\B,S 



with 



Lt .^ / II vec(xa;*)||B,sML(ds,da;), A := max(5R(cr(S))) 
Jo Jr-^ 



K2 B '■= max 



1X1 



xes+,\\x\\2=i\\\ycc{X)\\B,s , 
is the volatility process of a univariate MUC0GARCH(1, 1) process and y satisfies 

||vec(yt)l!i3,s<yt for allteR+ a.s. (4.2) 
Moreover, K^^b < ll^lll max^es+,||x||.=i(F5OT) ^ H^H^' 



88 



R. Stelzer 



Remark 4-2. (i) Provided S is unitary, K2,b = 1- Otherwise, an inspection of the proof 
shows that the inequahty (4.2) also holds if i^2,B is replaced by 115*112 in (4.1), which 
saves one from calculating the value of K2,b in practice. Likewise, \\A(E) A\\b,s can be 
replaced by \\A (g) = \\A\\l since ||(A (g) A){{C + Yr,-)^/^ ® (C + Yr,-)^/^)\\B.s < 
||5||i||5-i||i||(A^A)((C + yr,-)^/^®(C + yr,-)^/^)||2. 
This can be done in all of the forthcoming results involving K2.B or ||A®yl||s,s- 
(ii) As can be seen from the proof, the diagonalizability of B is essential and, unfor- 
tunately, it seems very intricate to extend the result to the non-diagonalizable case. In 
applications, however, this appears to be no severe constraint, as the non-diagonalizable 
matrices have Lebcsgue measure zero. 

Since the finiteness of moments of univariate C0GARCH(1, 1) processes is well known 
from [26], Section 4, we can now give sufficient conditions for the MUCOGARCH volatil- 
ity process to have some finite moments, which we will improve upon in Proposition 4.7. 

Proposition 4.3. Let k & Yq G § such that i^dlFoll'^) < 00 and let B he diagonaliz- 
able. Assume further that the MUCOGARCH volatility process Y is driven by a Levy pro- 
cess L satisfying E{\\Li\\'^'') <oo. Then E{\\Yt\\'') < 00 for allteM.+ and t^ E{\\Yt\\'') 
is locally bounded. 

4.2. Markovian properties and stationarity 

Turning to the study of the Markovian properties of a MUCOGARCH process, we refer to 
standard references like [15, 17, 18] for definitions and necessary general results. Moreover, 
we implicitly assume that our given filtered probability space is enlarged as in [36], 
page 293, to allow for arbitrary initial conditions of the SDEs, and the weak Feller 
property is defined as in [15], namely by demanding that the transition semigroup is 
stochastically continuous and maps the bounded continuous functions on the state space 
into themselves. 

The usual results on the Markov properties of SDEs (see [36], Section V.6) extend 
to locally Lipschitz SDEs on open sets (see [40], Section 6.7.1.2, for details) and to 
closed sets, provided the solution is ensured to stay in the closed set at all times and 
the SDE is defined on an open set containing the closed set. The latter is the case for 
the MUCOGARCH, the closed set and the open set Uc,e as defined in the proof of 
Theorem 3.2. Thus, we obtain the following result. 

Theorem 4.4. The MUCOGARCH process (G, Y) and its volatility process Y alone are 
temporally homogeneous strong Markov processes on R'^ x §J and Sj", respectively, and 
they have the weak Feller property. 

One of the most important questions regarding Markov processes in applications is the 
existence of stationary distributions. 
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Theorem 4.5. Let B e Afd(M) he diagonalizahle with S G GLd{C) such that S~^BS 
is diagonal. Furthermore, let L be a d-dimensional Levy process with non-zero Levy 
measure, X be defined as in Theorem 4-1 and a\ := ||S'||2||5'~^|||if2,B||^<8)^||_B,s- Assume 
that 

f log{l + ai\\vcc{yy*)\\B.s)Mdy)<-2X. (4.3) 

There then exists a stationary distribution fj, G A^i(§J), that is, the set of all probability 
measures on the Borel-a -algebra of Sj , for the MUC0GARCH(1, 1) volatility process Y 
such that 

f {{l + ai\\vcciyy*)\\B.st -'L)Mdy)<~'2Xk (4.4) 
Js.'' 

for some k eN implies that L+ ||a;||'^/i(dx) < oo, that is, that the kth moment of the 

d 

stationary distribution is finite. 

Of course, this result immediately translates to stationarity of V . Moreover, for c? = 1, 
it recovers the necessary and sufficient stationarity condition of [26]. 

Remark 4-6. (a) From [26], Lemma 4.1(d), it follows that, if (4.4) is satisfied for k e N, 
then it is also satisfied for all fc G N, k <k. 

(b) Combining the results of Section 4.1 shows that ai = H^Hi and || • \\b.,s = II ■ II2 if 
B is normal. 

Establishing uniqueness of the stationary distribution and convergence to the station- 
ary distribution for arbitrary starting values appears to be a rather intricate question 
due to the Lipschitz property holding only locally and the fact that d[L,L]'' lives on the 
rank one matrices. However, in the next section, we obtain at least asymptotic second- 
order stationarity and that the stationary second-order structure is unique under some 
technical conditions. 

To conclude this section, we consider some examples exploring the dependence of the 
stationary distribution on the parameters and the relation to the stationarity of univariate 
COGARCH processes. 

Example 4-1- Let cG M+\{0}, A,Be Mrf(M) and L be a d-dimensional Levy process. 
If V satisfies 

dVt = iB{Vt^ - cLd) + (Vt- - cId)B*) dt + AVll^ d[L, L\\ V^l^ A\ (4.5) 

then Z defined by Zt~Vt/c satisfies 

AZt = {B{Zt- - Id) + {Zt- - Id)B*)dt + AZ]I^ d[L, L\\ z'J^ A\ (4.6) 

which does not depend on c. In particular, if /i G A^i(§j") is a stationary distribution 
for (4.6), then /^^ G 7Wi(S+), defined by /ic(VK) = [liy^jc) for all Borel sets 1^ C Sj^, is a 
stationary distribution for (4.5). 
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Example 4-2. Assume that A,B,C and Yq are diagonal, the components of Yq are in- 
dependent and the components of the Levy process are completely independent, that is, 
whenever L has a jump, then only one of the d components jumps. In this case, Y or V, 
respectively, consists of d independent univariate C0GARCH(1,1) volatility processes. 
If each of the d univariate C0GARCII(1, 1) volatility processes converges in distribu- 
tion to a stationary distribution, then Y oi V, respectively, converges in distribution to 
a stationary distribution. In this example, condition (4.3) can be shown to imply the 
necessary and sufficient stationarity condition of [26], Theorem 3.1, for all components 
simultaneously. Actually, condition (4.3) is stronger than requiring that the univariate 
stationarity condition be satisfied for all components. 

However, it should be noted that the picture is very different when Yq is not diagonal 
because then jumps in one component of L typically affect all components of Y . Hence, 
it is not clear whether one still has convergence to a stationary distribution and whether 
this has to be the same distribution as the limit distribution when Yq is diagonal. When 
we have that Y is asymptotically second order stationary (see Theorem 4.20 below) and 
the limiting distribution for a diagonal Yq has finite second moments, the off-diagonal 
(covariance) elements of Y or V, respectively, necessarily converge to zero in as t ^ oo. 

Example 4-3- A degenerate situation occurs if C G Sj\§j"^. Take d = 2, A = al2, 
B = -ph with a,/3 e ]R+\{0} and C = J). If one has that Yq = yC with y e R+ 
(possibly random), then Y is at all times a scalar multiple of C and, when L jumps, all 
components (variance and covariance ones) have a jump of the same height. However, 
one again has a completely different picture if Yq is chosen differently, for example, Y is 
in S^"'" at all times, provided Yq G Sj"^. This also shows that the assumption C G Sj"^ 
made in the definition of the MUCOGARCH processes is essential to avoid pathological 
situations. 

4.3. Second-order moment structure 

Assuming stationarity and the existence of the relevant moments of the stationary solu- 
tion, we calculate explicit expressions for the moments of a stationary MUC0GARCH(1, 1) 
volatility process in this section, treat the non-stationary case along the way and present 
results regarding (asymptotic) second-order stationarity. Due to the special structure of 
the stochastic differential equation (3.4), especially due to the presence of the matrix 
square root, it is only possible under certain assumptions on the Levy process to obtain 
explicit formulae. The results of this chapter provide the basis for method of moments 
estimation, provided the volatility process is (approximately) observed and shows that 
the second-order structure of the volatility process is in line with observed financial data, 
since the matrix exponential decay of the autocovariance is rather flexible and has been 
found realistic in the analysis of OU-typc models (see [34]). 
Henceforth, we often assume the following in this section. 

Assumption 4-1 • (^t)teM+ is a second-order stationary yi\!iGOGAllGVi{\,\) volatility 
process. 
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Assumption 4-2. The pure jump part of the driving Levy process (it)tgR+ has finite 
variance which is a scalar multiple of the identity matrix: 

If we let Lf := J^^^^^^^x{fiL{ds,dx) - dsi^Lidx)) + J^^^^^^^xfiL{ds,dx) denote the 
pure jump part of L, then this means that there exists a ul ^ such that Yai{L\) = 
/jjd xx*VL{dx) ^(TLh- 

This assumption is comparable to considering only standard Brownian motion in 
Brownian-motion-based models and, hence, not very restrictive since any Levy process 
with finite second moments can be transformed into one satisfying Assumption 4.2 by a 
linear transformation and since the variance of G can still be flexibly modelled via the 
remaining parameters, as will be seen from Proposition 5.2. 

First, we need a refinement of Proposition 4.3 to the case where B is not diagonalizable. 

Proposition 4.7. LetY 6e a MUC0GARCH(1, 1) volatility process andk^{l}\j[2,oo) . 

If E{\\Yn\\^) < oo and E{\\Li\\^'') < oo, then E{\\Yt\\'') < oo for all teR+ and t ^ 
E{\\Yt\\'') is locally bounded. 

We can now calculate the expected value of the volatility. 

Theorem 4.8. Assume that Assumption 4-2 holds: 

(i) //£;(||ro||) <oo, then 

E{Yec{Yt))=e'^'E{vec{Yo)) + /* e-®(*~^Ms (7l(A A) vec(C) 

Jo 

with S§ := B (g) Id + Id® B + (TlA® A. If SS is invertible, then 

E{vec{Yt)) = e-®*(£:(vec(ro)) + ^^^-^(A ® A) vec(C)) - crL^-^{A ® A) vec(C) (4.7) 
for allteR+. 

(ii) Under Assumption 4-1, the stationary expected value E{Yq) of the MUCOGARCH 
volatility process satisfies 

BE{Yo) + E{Yo)B* + aLAE{Yo)A* = -ctlACA*. (4.8) 

// ^ is invertible, then the following formulae hold: 

E{vec{Yo)) = -aL.'^-\A ® A) vec(C) and 

(4.9) 

E{veciVo)) = ^-\B (g>Id + Id®B) vec(C). 

Remark 4-9. Observe that the stationary expectation is the limit of the expected value 
in (i) for t oo provided ct(.^) C (-oo, 0) + iR. 
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Theorem 4.5 can not only be used to show that Assumption 4.1 is satisfied, but also 
to ensure the invertibility of Hence, Theorem 4.8 provides an explicit expression for 
the mean of the stationary distribution of Theorem 4.5. 

Lemma 4.10. Assume that (4-4) "is satisfied with k = 1 for the MUCOGARCH volatility 
process Y and that Assumption 4-2 holds. Then, SS, defined as above, is invertible and 
cr(^) C (-00,0) + iK. 

Before analyzing the variance, let us study the autocovariancc function. If {Xt)t£R+ 
is a second-order stationary process with values in M'', the autocovariancc func- 
tion acovx:K Md{M.) of X is given by acovx(/i) = cov{Xh,Xo) = E{XhXQ) - 
E{Xo)E{Xo)* for h>0 and by acovx(/i) = (acovx(— /i))* for h < 0. As we arc con- 
sidering matrix-valued processes {Zt)teM in the following, we set acov^ := acovvcc(z) in 
this case. 

Theorem 4.11. (i) Under Assumptions 4-1 and 4-2, the autocovariance Junction of the 
MUCOGARCH volatility process satisfies 



acovy(/i) = acovv(/i)=e(-^®-f''+^''®-^+'"^^®-^)''var(vec(yo)), h>0. (4.11) 



(ii) If Assumption 4-2 is satisfied and £'(||yo||^),i?(||ii||^) are finite, it holds that 



cov(i;+,,y„) = cov(K+/., K) = e(^»^-+^-®^+'^-^®^)'' var(vcc(r„)) (4.12) 



for all u,h>0. 

The autocovariance function of the volatility process Y is thus exponentially decreas- 
ing in a matrix sense, so the individual entries may decay as sums of exponentials, 
exponentially damped sinusoids (if the eigenvalues have non- vanishing complex parts) or 
exponentially damped polynomials (if the matrix is not diagonalizablc) . 

However, we are so far lacking an explicit expression for var(vec(yo)). Unfortunately, 
our Assumption 4.2 on the second moment of the jumps of the driving Levy process L 
seems not to be sufficient to obtain an explicit expression for the variance. 

As we shall see from the proofs, the quadratic variation of the vectorized discontinuous 
part of the quadratic variation of the driving Levy process. 



— acovy (/i) ^{B®Id + Id®B + glA ® A) acovy (/i) 



(4.10) 



forh>0. 
Hence, 



[vec([L,L]° 



),vcc{[L,L]^)]^ — / / vcc(a;a;*) vec(a;a;*)*^_L(ds,da;), 



JO JR'' 
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which is again a pure jump Levy process of finite variation, will appear in our calculations 
of the second moment and we need it to have finite expectation. In fact, we even need to 
make specific assumptions on its expectation. 

To determine what assumptions are reasonable, let us assume for a moment that L 
is a d-dimensional compound Poisson process with rate one and with jump distribution 
being the c?-dimensional standard normal distribution. This implies that [-Zj,-/^]'' is a 
compound Poisson process with rate one and with the jump distribution being a Wishart 
distribution. Then, denoting the d-dimcnsional standard normal distribution by TV (da;) 
and noting that vec(a;a;*) vec(a;a;*)* = {x ^ x){x* (8) a;*) = (xx*) (8) {xx*), we have 



from [30], Theorem 4.1. Here, Kd S Md2{M.) denotes the commutation matrix which can 
be characterized by Kdvec{A) = vec{A*) for all A E Md{M.) (sec [30] for more details). 
This generalizes to the following result. 

Lemma 4.12. Let L be a d- dimensional compound Poisson process with rate c and 
with jumps distributed like \fsX , where X is a d- dimensional standard normal random 
variable and e is a random variable in with finite variance and independent of X . 



E{[vec{[L,L]^),vec{[L,L]y]l) = cE{e^){Id2 + Kd + vec{Id)yec{Idy)- (4-14) 



Moving away from a Levy process of finite activity, a similar result holds for the 
following variant of type G processes, a special kind of a normal mixture. 

Definition 4-13 (Type G). Let L be a d-dimensional Levy process. If there exists an 
W'^ -valued infinitely divisible random variable e independent of a d-dimensional standard 
normal random variable X such that Li = y/eX , then L is said to be of type G. (Here, 
= denotes equality in law.) 

We have chosen the name "type G" above because these processes correspond to a 
particular case of multG laws as defined in [2], Definition 3.1. Actually, many interesting 
Levy processes are of type G, for instance, the multivariate symmetric GH (NIG) pro- 
cesses with the parameter E set to Id (see [9, 31]). For details on distributions and Levy 
processes of type G in general, we refer to [2, 29]. 



Lemma 4.14. Let L be a d-dimensional Levy process of type G with a finite fourth 
moment. Then E{[vec{[L, L]^),vec{[L, L]°)]l) = pL^Id^ + Kd + yec{Id)vec{Id)*) with pL € 




(4.13) 



/rf2 + Kd + vec(/rf) vec(/d)* , 



Then, 



+ 



These results motivate the following assumption. 
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Assumption 4-3. The pure jump part of the driving Levy process (it)tgR+ has a finite 
fourth moment, that is, /jj^ ||x|j*i^L(da::) < oo, and there exists a real constant such 
that 

S([vec([L, L]^), vcc([L, L]^)]?) = pl(/<j2 + + vcc(/<j) vcc(/d)*). 

Intuitively, this means that the jumps of L have the same fourth moment as a standard 
normal distribution. This assumption is considerably more restrictive than Assumption 
4.2. However, it is comparable to the assumptions made for discrete-time multivariate 
GARCH processes (see [21]) and from the proofs, one sees that exphcit results are only 
obtainable if the fourth moment of the jumps is comprised of well-understood matrices 
which act on tensor products in a suitable way. 

To state our next result, we need to introduce some additional special linear operators 
and matrices. If we define 

Q:Md2(R)^Md2(R), 

(4.15) 

{QX){k-l)d+L{p~l)d+q — X{k-l)d+p,{l-l)d+q for aU fc, l,p, g — {1, 2, . . .,d}, 

then = Q, obviously, and Q(vec(X) vec(Z)*) = AT (g) Z for aU X,Z &Sd (see [34], 

Theorem 4.3). Furthermore, we define Q € M^4(R) as the matrix associated with the 

linear map vccoQ o vcc^ on M and JCd S Mrf4(IR) as the matrix associated with the 

linear map vec{Kd vec~^ (x)) for x eM.'^ , where vec : M^^ (R) — R'^ .It is easy to see that 

both Q and K-d simply permute the entries of a vector x S M . Both Q and K-d are thus 

permutation matrices, so we have ||Q||2 = ||^d||2 = 1, where || • II2 is the operator norm 

,4 

associated with the usual Euclidean norm on R . 

Theorem 4.15. Assume that Assumptions ^.2 and 4-3 hold. 
(i) //£;(|lyo|l^) <oo, then 

^vec(£;(vec(rt)vec(yt)*)) 
^^Eivcc{Yt)®vec{Yt)) 

(4.16) 

= <rvec(£;(vec(rt) vec(yt) )) + (ctl(A® A) ® 7^2 + £/n)vec{C) ® £;(vec(yt)) 
+ (crL/rf2 (g) (AgjA) +i^/7^)^;(vcc(Yt)) g) vec(C) +£/7^vcc(C) (g) vec(C), 

where 

£/ =(A(gA)(g(A(gA); n = pL{Q + ICdQ + Id^) 
'rf {B <g) Id + Id ® B) <g) 7^2 -|- 7^2 ®{B®Id + Id®B) 
+ aL{{A® A)® Id2 +Id2(»{A(E)A))+ ^/7^. 
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(ii) Under Assumption 4-1, the stationary second moment i?(vec(lo) vec(lo)*) of the 
MUCOGARCH volatility process satisfies 

.^E{YCc{Yo)vcc{Yoy) + E{vcc{Yq)vcc{Yo)*)^* 
+ {A(g) A)RE{vcc(Yo)YCc{Ya)*){A* (g) A*) 

(4.17) 

= -(Jl[{A (E) A) vcc{C)E{vcc{Yo))* + E{vcc{Yo)) vcc{C)*{A* A*)] 

~{A(E) A)Ii{E{vec{Ya)) vec(C)* + vec(C)£;(vec(ro))* + vec(C) vec{C)*){A* ® A*), 

with R := Pi (Q + KdQ + Id^)- 

Provided is invertihle, E{ycc{Yi))ycc{Yq)*) is given by 

vec(£;(vec(ro)vcc(yo)*)) 

= -'^#'-^[£/7^(vec(C) ® vcc(C)) 

(4.18) 

+ (o-L {A<SiA)® 1^2 + .s/n) vec(C) ® £'(vcc(ro)) 
+ {(TLld^ «) {A(gi A) + .s/n)E(vec{Ya)) (gi vec{C)]. 

Remark 4-16. The differential equation (4.16) is an inhomogeneous linear differential 
equation with constant coefheients. Hence, it is standard to obtain an explicit solution. 
We refrain from stating it, as the stationary case seems to be of the most importance. 

Again, condition (4.4) of Theorem 4.5, which ensures the existence of moments of the 
stationary distribution obtained there, also implies invertibility of under an additional 
technical assumption. 

To state the result, we set S = Sg)S®S^S and define a new norm || • ||^^ on 

M'' by setting H^^H = ||'5^"'^a;||2. The associated operator norm on Afrf4(M) is given by 

\\x\\^g = \\s-'xsh. 

Lemma 4.17. Assume that (4-4) satisfied with k ^2 for the MUCOGARCH volatility 
process Y and that Assumptions 4-2 o.'i^d 4-3 hold. Provided that 

\\Q + }CdQ + W^s ^ KlB\\Yec{Id2 +Kd + vec{Id) vec(/d)*)|l^ (4.19) 

also holds, o{f;€) C (— oo,0) + iM and is invertihle. 

A rather unpleasant feature of this lemma is that we need the technical condition 
(4.19). The following lemma shows that it is always true if S is unitary, and for concrete 
parameter values, it can, of course, be checked numerically. 

Lemma 4.18. If S is unitary, then (4-19) is satisfied. 

We end our comprehensive calculations for the second-order moment structure of the 
MUCOGARCH volatility process by turning to the stationary variance. 
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Corollary 4.19. If Assumptions 4-l~4-3 hold and 'if is invertihle, then the stationary 
variance var(vec(yo)) = var(vec(Vo)) of the MUCOGARCH volatility process is given by 

vec(var(vec(yo))) = -'^'^[{crl'^i^-^ <E) + s^TZ){vec{C) (E) vec{C)) 

+ {aL{A^ A) ® 1^2 + £/n)vec{C) ® E{vec{Yo)) (4.20) 
+ (f7L/d2 ®{A®A) + £/TZ)E{vcc{Yo)) ® vcc(C)] . 

Proof. Combine (4.18), 

vec(i;(vec(yo))£^(vcc(yo)*)) = crK^"^ i^-i)^/(vcc(C) vec(C)) 

and the elementary formula var(vec(Fo)) = E{vcc{Yo)vcc{Yo)*)- E{vcc{Yo))E{vcc{Yo)*). □ 

Under specific moment assumptions on the driving Levy process, we have thus calcu- 
lated the second-order structure of a stationary MUCOGARCH volatility process com- 
pletely. 

Finally, we give conditions ensuring (asymptotic) second-order stationarity. A stochas- 
tic process X in is said to be asymptotically second-order stationary with mean 
/i e M'* , variance E £ §^2 and autocovariance function / : K+ Mj2 (M) if it has finite 
second moments and 

lim E{Xt) = fJ,, lim var(vec(Xt)) = S and 

lim sup {II cov(vec(Xt+,,),vec(Xt)) - /(/i)||} = 0. 

Theorem 4.20. Let Assumptions ^.2 and 4-3 be satisfied and assume further that the 
matrices B,3§,''^ are such that a{B),a{3§),a{'^) d (-00, 0) 4- iR. 

(i) If Yq satisfies (4.9) and (4.20), then the MUCOGARCH volatility process Y is 
second- order stationary. 

(ii) // i?(||yo|P) < 00, then the MUCOGARCH volatility process Y is asymptotically 
second-order stationary with mean, variance and autocovariance function given by (4-9), 
(4.20) and (4.11). 

5. The increments of the MUCOGARCH(l, 1) 
process 

Thus far, we have mainly studied the MUCOGARCH volatility processes Y and V . 
However, in practice, one typically cannot observe the volatility, but only the process G 
(which, in a financial context, for instance, resembles log-prices) at finitely many points in 
time. In the following, we presume that G is observed on a discrete-time grid starting at 
zero and with fixed grid size A > 0. It is obvious how the results of this section generalize 
to non-equidistant observations or to the setup considered in [10, 26]. 
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In financial time series, one commonly observes that the returns themselves are un- 
correlated, but the "squared returns" (that is, the return vector times its transpose in 
a multivariate setting) are considerably correlated. The following results show that the 
MUCOGARCH model can reproduce this very important stylized feature and, further- 
more, they provide the basis for simple moment estimators (as in [22]). 

We define the sequence of increments G = (G„)„gN by setting 

G„= / K-'dL,. (5.1) 

J(n-1)A 

Moreover, we shall throughout most of this section presume the following. 
Assumption 5.1. Y (or equivalently V ) is stationary. 
Proposition 5.1. If Assumption 5.1 holds, then G is stationary. 

Proof. Employing Theorem 4.4 and the same arguments as for [26], Corollary 3.1, shows 
that G has stationary increments. □ 

In order to be able to obtain explicit expressions for the moments of G, we need to 
strengthen Assumption 4.2 as follows. 

Assumption 5.2. Assumption ^.2 is satisfied and, moreover, 

E{Li) — and var(Li) — {aw + cfL)Id with uw > 0. 

This assumption means that in addition to Assumption 4.2, the Brownian part of L is 
a scalar multiple of d-dimensional standard Brownian motion. 

We start by giving conditions for the finiteness of the second moments of G, and thus 
of G, without requiring stationarity and explicit expressions for the moments in the 
stationary case. 

Proposition 5.2. Assume that E{Li) = 0, E{\\Li\\'^) < oo and E{\\Yo\\) < oo. Then 
E{\\Gtf) <oo for allteR+. 

If Assumptions 5.1 and 5.2 are also satisfied, then the stationary sequence G has the 



following second-order structure: 

S(Gi)=0, (5.2) 

var(Gi) = [aL + aw)AE{Vo), (5.3) 

vec(var(Gi)) = {aL+aw)A^~HB(^Id + Id(^B)vec(C), (5.4) 

acovG(/i) = for all h e Z\{0}. (5.5) 
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Remark 5.3. (i) In the stationary case, this shows that G is a white noise sequence. 

(ii) Straightforward extensions of the arguments in the proof show that if Y is not 
stationary, but only (asymptoticaUy) second-order stationary, then G is (asymptoticahy) 
second-order stationary. 

For the squared returns GG* ~ (G„G* )„gN, we get the following. 

Proposition 5.4. Assume that E{Li)=Q, S(||Li||'^) < oo and EiWYoW^) is finite. Then 
E{\\Gt\\^)<(x, and, likewise, E {\\G tGlW^) <qo for allteR+. 

If Assumptions 5.1 and 5.2 are also satisfied, then the stationary sequence GG* has 
the following second-order structure: 

E{GiGl) = {aL+<yw)^E{Vo), (5.6) 
acovGG-(/i) -=c^^''^-i(/<j2 -c-^^)(crL + fTvi/)cov(vec(yA),vec(GiG*)) (5.7) 

for heN. 

Thus, the squared returns GG* have, like an ARMA(1, 1) process, a matrix exponen- 
tially decreasing autocovariance function from lag one onwards. That such an autoco- 
variance structure is reasonable for financial data can be seen from [34], for instance. In 
financial data, (quasi-) long-range dependence is frequently encountered, which is often 
(see, for example, [4]) well modelled by specifying the autocovariance function of the 
squared increments as the sum of fast and very slowly decaying exponential functions. 
Since we have a matrix exponential decay, we obtain such a behavior componentwise 
by appropriate choices for our parameters, with the different rates of the exponential 
decay being determined by the eigenvalues of Additionally we can cover a sinusoidal 
component. 

In the univariate case, [22] obtained, under additional assumptions on L, explicit ex- 
pressions for var(vec(GiGJ)) and cov(vcc(Ya), vec(GiGJ)). As these arc, however, al- 
ready rather lengthy and complicated formulae, we refrain from calculating these values 
in our multivariate model. 

6. Proofs and auxiliary results 

In this section, we provide the proofs of our results, along with necessary additional 
technical results. 

6.1. Proofs for Section 3 

We begin with some matrix analytic results analyzing the Lipschitz properties of the map 
V ^ (g) used in the definition of the MUC0GARCH(1, 1) volatility process. We 
denote by || • [[2 the operator norm associated with the usual Euclidean norm on R''. 
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Lemma 6.1 ([7], Problem 1.6.11). For all A, Be Afd(R), we have 
\\A(g, A~ B(g) B\\2<2ma.x{\\A\\2,\\B\\2}\\A- B\\2. 

In particular, the mapping (g) : Md(R) — > M^2(M.), X X (I) X is uniformly Lipschitz on 
any set of the form {x G A/d(R) : ||a:|| < c} with c > 0. 

The proof is obvious from the ideas outhned in [7] . 

Lemma 6.2 ([7], page 305). Let e §J and a > such that A,B> aid- Then 

\\A'/'-B'/%<-^\\A-B\\2. 

Hence, the mapping — > i— )■ X^^"^ is uniformly Lipschitz on any set of the form 
{xe§j:x> cid} C S,j+ with c> 0. 

For a variant of the above statement see [24], page 557. 

Lemma 6.3. Consider the map F:§+ ^ §+,X i-^ X^/^ ^ ^1/2 ^ (x ® X)^/'^. F is 
continuous and uniformly Lipschitz on any set of the form {x G E>'^ : x > cl, \\x\\ < c} with 
c,c> 0. Moreover, we have that WA^/"^ <Si A-^^^IU = \\A\\2 for all A€§d- 

Proof. The identity X^/^ ® X^/^ = (X (g X)^^^ is an immediate consequence of basic 
properties of the tensor product (see [24], Chapter 4) and the continuity of F follows 
from the continuity of the tensor product and the positive definite square root (see [24] , 
Theorem 6.2.37). The Lipschitz property follows from a combination of the previous two 
lemmas. Finally, \\A^/'^ (g) A^^'^y = \\A\\2 is established by noting that || A^/^ ^ ^i/Zjj^ = 
(see [7], page 15) and \\A^^^\\2 = \\A\\l^^ . The latter follows immediately from 
the fact that \\A\\2 = ^p{A*A) = p{A) for all ^ e §+. □ 

Finally, we show that the global Lipschitz property is not satisfied for this map, not 
even if we restrict it to sets being bounded away from zero. 

Lemma 6.4. For the map F defined in the previous lemma, there exists no finite K e 
such that 

\\F{x)-F{y)\\<K\\x-y\\ (6.1) 
for all x,y £ Sj^ . The same holds for all x,y £ {z eSd'- z > C} with arbitrary C G . 

Proof. From the following proof, it is clear that we can take d ~ 2 without loss of 
generality. Let x = diag(a;i, 3:2) and y = diag(yi, 2/2), with xi, 0:2,2/1, 2/2 G R^\{0} and 
diag(xi,X2) being, as usual, the diagonal matrix with diagonal entries xi and X2- We 
have F{x) = diag(a;i, ^a;iX2, ^/xiX2, X2)- Assume that (6.1) is true with a finite K G K"*". 
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There is then a finite k G R+ such that \y^xiX2 — ^/WlBl < k{\xi — yi \ + \x2 ~ J/2I) for 
aU a;i,a;2,2/i,2/2 e M+\{0}. Choosing X2 = y2 = 1, this gives \^Jxl - y/yi\ < - yi\ 
for all Xi,yi G R+\{0}, which is a contradiction to the well-known fact that the square 
root is not globally Lipschitz on R+\{0}. Regarding the case x, j/ G {z G Sd : z > C}, we 
can, without loss of generality, restrict ourselves to C = cid with c G R"'"\{0}. Choosing 
2^2 = 2/2, xi~ 9c and yi = 4c gives |-y/cx2| < hkc. As X2 can be taken arbitrarily large. 



In the following, we use the fact that any stochastic differential equation defined on an 
open set which has locally Lipschitz coefficients growing at most linearly has a unique 
solution until the first time the open set is left or its boundary is reached. This result 
follows along the same lines as the usual existence results for SDEs with locally Lipschitz 
coefficients defined on (see, for example, [33] or [36], Theorem V.38). Alternatively, 
a proof for open sets of the type relevant below can be found in [40], Section 6.7, which 
uses only the standard existence and uniqueness results for SDEs on W'- with globally 
Lipschitz coefficients and orthogonal projections. 



Proof of Theorem 3.2. Define the maps F and G by F(vcc(y)) — {Id ® B + B ® 
/d) vec(j/) and G{y) ^{A® A){{C + yf^ ® (C + yf^). The SDE (3.4) can then be 



Moreover, we define the set Uc,e ~ {x <E'Bd'-x> —eld] for some e with < e < mincr(C). 
Then the set Uc.e (and thus vcc{Uc.e)) is open and for each x G Uc.e^ we have x + C > 
(mincr(C) — e)Id G ^. Since the foregoing results imply that G is locally Lipschitz on 
Uc,e and has linear growth (a function / has linear growth if ||/(2;)||^ < C(l + ||a;|p)), 
standard results on the existence of solutions of SDEs give that (6.2) has a unique locally 
bounded solution {Yt)t£s.+ with initial value loj provided it can be ensured that every 
solution does not leave the set Uc,e or touch its boundary. However, it is easy to see 
that every solution must satisfy Yt > c^^YqC^'* since all jumps AV^/^ A[L, L]^Vf^/'^ A* 
are positive semidefinite and between jumps, Y follows the deterministic differential 
equation dYt = (BYt- + Yt-B*)dt uniquely solved by Yt = c^*Yoe^'\ so any solution 
necessarily stays in Sj". 

The finite variation property is clear since time t and [L,L]'' are of finite variation. □ 



Proof of Theorem 3.6. Define Mt = A{C + Y,,-)^/^ d[L, L]l (C + Y,^)^^^A*. Then 



M is Sj'-incrcasing and of finite variation and Y obviously solves the stochastic differ- 
ential equation dXt = {BXt- -\- Xt-B*)At + dMt{*). Standard theory implies that this 
differential equation has a unique solution, and the same elementary calculations as for 
Ornstein-Uhlenbeck processes show that the solution of (*) with initial value Fq: which 
is necessarily equal to 1", is given by 



this is a contradiction. 



□ 



written as 



dvcc(rt) =F(vec(rt_))dt + G(yt-)dvcc([L,i]t^). 



(6.2) 
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= e^*roe^**+ re^(*-^)A(C + y,_)i/2d[i,i]°(C + y,_)i/2A*e^*(*-^). ^ 
Jo 



6.2. Proofs for Section 4 

6.2.1. Proofs for Section ^.1 

The univariate C0GARCH(1, 1) bounds wih first be shown for processes driven by com- 
pound Poisson processes and then for the general case using an approximation by com- 
pound Poisson processes which is of interest in its own right, as it provides, for instance, 
a possible approximation scheme to be used in simulations. 

To see that the processes defined in the following are indeed univariate COGARCH 
processes, we need the following general lemma. 

Lemma 6.5. Let {Lt)t(z^+ he a driftless Levy subordinator. There then exists a Levy 
proeess {Lt)t^TSi+ in R such that Lt ~ [L, L]^ for all t G R+. 

Proof. Denote the jump measure associated with Lhy hl, that is, Lt — /g a;/iL(ds, da;), 
and denote its Levy measure by vl- Let 

Lt= / Va;(A'L(ds,da:) - dsi'L(da;)) + / / V^^L(ds, da;), 

Jo Jo<x<l Jo Jx>l 

noting that the existence of the first integral follows from the finiteness of /q^^.<2 sf^ (da;) = 
Jp^^^j^ a;;/L(da;) since L is of finite variation. Now (Lt)tes& is a- Levy process and 

The following elementary results regarding the norm j| • ||b,s are straightforward to 
obtain. 

Lemma 6.6. It holds that \\S®S\\b,s = \\S\\l and \\S^^®S^^\\b.s = \\S^^\\l. Moreover, 

||a;||B,s<||5"'!l2l|a^ll2 and ||.Ti|2 < ||5||2||a;||B,s for all x ^W^\ 
\\X\\B,s<\\S\\l\\S-^\\l\\Xh and \\Xh <\\S\\1\\S-Y^\\X\\b,s for all X G Md^iR). 

Note that a very similar norm has also been used in [10]. 

We can now analyze the norms of compound Poisson driven MUCOGARCH volatility 
processes. Recall that in the univariate case, the MUCOGARCH volatility process Y is 
just a deterministically scaled version of the COGARCH volatility process Y defined in 
[10]. 



Proof of Theorem 4.1 if L is compound Poisson. Using Lemma 6.5, it is clear that 
the process {yt)tm defined in Theorem 4.1 is a univariate MUC0GARCH(1, 1) process. 
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Let Ti be the time of the first jump ofLand let f € [CTi). Since \\e^^-'<»B+Bmic)ty^^^-_ 
e^-^*, it holds that 

II vec(rt)||B.5 - ||c(^-«^+^®^'')* vec(ro)|lB s < Ije^^'^^^+^^^'^'i^ vec(yo)||B,s 



Thus, (4.2) is shown for aU t e [0,ri). At time Ti, we have 
l|vec(rrJb.5 

= llvcc(yri-) + {A<E>A){{C + Yr,-)^^^ <E>{C + Fn-)'/') vec(ALr, {ALrjnWB.s 

<yr,- + \\A(g> A\\BM\iC + Yr,-y/^ <E> {C + Yr,-y/^BM\^(=c{ALrA^Lrjn\\B,s 
<yr,- + \\A(^A\\BM\S\\l\\S-X\\{C + Yr,-)^/^(^{C + Yr,^y^''hALr, 
< yr,- + \\A®A\\B.s\\S\\l\\S-Y2{\\Ch + \\Yr,-h)ALr, 

<yri- + P®A||B,5||5||^||5"'||^i^2,B(/^2:sll^^ll2 + l|vcc(rr,_)||B,s)AZr,=yri, 

which establishes (4.2) for t = Fi. Iterating these arguments shows (4.2) for all t G M+. 

The first inequality for K2.B follows from Lemma 6.6 and the second one by [23], page 
314. □ 



In order to extend Theorem 4.1 to MUCOGARCH processes driven by general Levy 
processes, we need to show that we can approximate a MUCOGARCH volatility process 
by approximating the driving Levy process. The following result is very similar to [10], 
Lemma 8.2. However, we need to give a detailed proof since the standard results cannot 
be applied due to the fact that we have only locally Lipschitz coefficients. 

Proposition 6.7. Let Y be a MUCOGARCH volatility process with C £ §J+ and Yq g 

Sj, driven by a Levy process L in , and let (e„)ngN be a monotonically decreasing 
sequence in ]R+\{0} with lim„_>ooen =0. Define compound Poisson Levy processes Ln 
by Ln.t = /o /ijd ||x||>£ x^Li<is,dx) for n eN and associated MUCOGARCH volatility 
processes Yn by 

dy„,t = (By„,t_ + r„,t_B*) at + A{C + Fn.*.)^/' d[L„,L„]? {C + Y^^t.fl^A*, 

(6.3) 

Ynfi — lo- 

Then Yn ^Y as n —> 00 almost surely uniformly on compacts. 

Proof. First, observe that [L„, i„]f = /o/||2;||>£ a;a;*/iL(ds, dx) implies that [i„, L„]'' 
[LjL]'', as n — )■ 00 a.s. uniformly on compacts, and that [L„,L„]'' is monotonically in- 
creasing in n. 
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Since all processes involved are of finite variation, we can prove the claim with a 
pathwise approach. So, fix a; £ O and thereby one path. Let T € K"*" be arbitrary. The 
Gronwall inequality (see [36], Exercise 15, page 358) shows that 

\\YnA2<(^\\Yoh + \\A\\l\\Ch£ Ij^^^^^ 

< (llYoh + WAWlWChJ^ I ||a;|l2^i(ds,dx))ell^ll'^"-^«''ll-ll^'^-('i^^'i-)+2ll^ll=^, 

\\Yth<(\\Yoh + \\A\\l\\Chl^'' JjM^^^^ 

for all t e [0,r]. Since Yt > e^*Yoe^'\ Y„^t > c^'FoC^** and Yq is positive semidefinite, 
C + Y and (C + y,i)„gN all remain in one common compact set in S^"*" on [0,T]. Thus, 
when considering (3.4) and (6.3), we can regard the coefficients of these SDEs as being 
globally Lipschitz with a common Lipschitz coefficient. Thus, [36], Corollary, page 261 
after Theorem v. 11, implies that Yn{uj) — )■ Y{uj) uniformly on [0,T]. Note that, formally, 
the result of [36] is applied on the probability space given by the set {w}, the trivial 
(T-algebra {{a;},0} (which also gives the filtration) and the Dirac measure with respect 
to UJ. 

Since Ld Eil and T E IR+ were arbitrary, this completes the proof. □ 

Proof of Theorem 4.1 for general L. Let (F„)„gN be the sequence of compound Pois- 
son driven MUC0GARCH(1, 1) processes converging a.s. on compacts to Y constructed 
in the last proposition. For n e N, denote by 2/„ the univariate MUC0GARCH(1, 1) pro- 
cesses with ]| yec{Yn^t)\\B.s ^ yn,t for all t G M+. Then y„^t + i^j'^UC'll is a univariate 
C0GARCH(1, 1) volatility process as defined in [26], where it is denoted by af_^_. Since 
we only add more jumps in [L„, in]" when we increase n, it is straightforward to see from 
equations (3.3) and (3.4) in [26] that Hn+i.t > Vn.t for all 7i,Z e N and t € M"*". Moreover, 
defining the process y by 

dyt = 2Xyt- dt+\\S\\l\\S-^\\lK2.B\\A<E>A\\B,s(^^ + Vt-) dL*, 



K- 



2,B 



(6.4) 

yo = II vec(yo)liB,5, 

with Lt := /q II vec(xa;*)||B,s/^L(ds, da;), the same argument implies that yn,t < yt for 
aU n G N and t € R"*". Note that (it)tgR+ is a well-defined Levy process, because there is 
sl K > such that 

i\\vec{xx*)\\B,sAl),^Lidx)<K [ {\\xx*\\2Al)i^Lidx) ^ K [ (||.t||2 A l)i.L(dx) < oo. 

Passing to the limit n — )■ oo in || vcc(y„.t)||B.5 < yn.t < yt establishes || vcc(yt)||B^s < yt 
foralHGM+. ' □ 
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Proof of Proposition 4.3. Let y be the process constructed in Theorem 4.1. It then 
suffices to show that E{yf) < oo and that this is locahy bounded in t. By construc- 
tion, E^yg) is finite. Moreover, let L be the Levy process constructed in Lemma 6.5 
such that Lt = [L,L]f. The finiteness of i?(|lLi|p'') implies that /jj^ ||a;||2''j^L(da;) = 
/jjd ||xx*|l2VL(dx) < oo (with || • ||2 denoting the Euclidean norm in the first integral and 
the associated operator norm in the second integral) . Since the finiteness of the integrals 
is independent of the particular norm used, it follows that J^^ \\vec{xx*)\\'^ ^^^{dx) = 
J^\x\''i'i{dx) = /jg |xp'^z/£(dx) < oo. Hence, E{\Li\'^'') is finite and using the results of 
[26], Section 4, as in the proof of [10], Proposition 4.1, completes the proof of this propo- 
sition. □ 



6.2.2. Proofs for Section 4.2 

In order to show the existence of a stationary distribution of the MUCOGARCH volatility 
process Y, we need to recall a result from the theory of weak convergence. For more details 
and the relevant background, we refer to any of the standard texts (for example, [8, 25]). 
Below, we denote by Aii{E) the set of all probability measures on the Borel a-algebra 
of a Polish space E. 

The following theorem on the existence of a stationary distribution for a Markov pro- 
cess is referred to as the "Krylov-Bogoliubov existence theorem" in the literature. For a 
proof, sec [13], Section 3.1, or [37], Theorem 4.6. 

Theorem 6.8. Let E be a Polish space and (Ps)sga+ the transition semigroup of an 
E -valued weak Feller Markov process. Assume that there is an rj € Mi{E) such that the 
set {Pf*rj:t £ R+} is tight. There then exists a ^j, € A4i{E) such that P^n = jj, for all 
t G K^, that is, ^ is an invariant measure for {Ps)seR+ or a stationary distribution for 
the Markov process, respectively, and ^ is in the closed ( with respect to weak convergence ) 
convex hull of {PtV ■ t € M^} . 

Above, Pt* -.MiiE) -J> Mi{E) denotes the operator given by Pt*^(J7) = /g Ps{x, U)^{dx) 
for any Borel set U where Pt{x,U) is the transition probability of the Markov process 
from the initial state x to the set U at time t G M+ . 

Proof. Proof of Theorem 4.5 Let X,L be defined as in Theorem 4.1 and L be the Levy 
process constructed in Lemma 6.5 such that Lt = [L,L]^. Then 

/ log(l + ai||vcc(yy*)||i3.s)i'L(dy)= / log(H- aiy^)i^£(dy) 

and thus [10], Theorem 3.1, (sec also [26], Theorem 3.1), shows that the process y satis- 
fying (4.1) converges in distribution to a distribution concentrated on R"*". Assume now 
that yo has this stationary probability distribution and is independent of (Ls)ggjj+. Set- 
ting Yq = II vcc(/rf)| | i; s ^'^ gives an initial value for the MUCOGARCH volatility process 
that is independent of L and, moreover, || vcc(yo)||B.s = 2/o- Thus, the process y satis- 
fying II vec(lt)||B_s < yt for all t G (see Theorem 4.1) is stationary. Since for every 
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K>0, the set {.t e Sd : |ja;|j < K} is compact in §+, P{\\Yt\\B,s <K)> Pijjt < K) and y 
is stationary with a stationary distribution concentrated on R+, it foUows that the set 
{^{Yt):t€ M+} of laws if (y*) of Yt forms a tight subset of Mi{^^). Therefore, The- 
orem 6.8 combined with Theorem 4.4 imphes that there exists a stationary distribution 
A^i(Sj) for the MUCOGARCH volatihty process Y such that ^ is in the closed 
convex hull of {if (Yt) : t G M+} . 

If (4.4) holds for some /c S N, [10], Proposition 4.1 (see also [26], Section 4), shows 
that the stationary distribution of y has a finite fcth moment. This, in turn, implies that 
-E'dl^ijl'^) — ^ fo'" some finite c S K+ and all t G M+. Hence, /g+ ||x||''"/i(da;) < oo because 
fj, is in the closed convex hull of {if (Yt) : t G M+}. □ 

6.2.3. Proofs for Section 4.3 

In the following calculations of moments of F, we often use the fact that the stochastic 
continuity of L implies E{Yt-) = E{Yt), var(vec(yt_)) = var(vec(yt)) and similar results. 
Moreover, the following version of the so-called compensation formula is needed. 

Lemma 6.9. Assume that {Xt)ti^^+ is an adapted cadlag M^C^) -valued process sat- 
isfying E(\\Xt\\) < oo for all t € R^, 1 1-^ E{\\Xt\\) is locally bounded and (it)tgK+ is a 
driftless pure jump Levy process in R'^ of finite variation with finite expectation E'dlLiU). 
Then E{jl X^- dL,,) = J* E{X,^)E{Li) ds for teR+. 

Proof. Since Lt = /J J^^ z/^i(ds,dz), the compensation formula (see [28], Section 4.3.2) 
implies that 

e( f X.^dL,) ^e( f f X,_ztiLids,dz)]^E( f X,^ f zz/i(dz)ds 

= e(^J^ Xs-E{Li)ds 

Observing that /j£;(|lX,_||)|l£;(Li)|l ds is finite for every t G M+ , an application of Fu- 
bini's theorem completes the proof. □ 

Proof of Proposition 4.7. Consider the case fc = 1 first. Elementary arguments give 
that 



\Yth< WYoh + W A\\i\\C\ 



f f ||a:||2^L(ds,dx) + 2i|Bi|2 / ||n-||2ds 
Jo Jm'' Jo 



\A\\l f [ i|r,„||2|lx-||^/.L(ds,dx-). 

JO JR'' 
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Using stochastic continuity, the compensation formula and the observation that 
/j[jd||a;||2'^L(da:)< oo due to the finiteness of i^dliiHj), this imphes that 



Emh) < E{\\Yoh) + WmChT / MiMdx) 



for aU te [0, T] and any T e IR+. The Gronwah lemma thus shows that i?(||yt||2) is finite 
and bounded for t E [0,T]. Since T was arbitrary, this completes the proof for this case. 

In the case /c > 2, we obtain from [36], Theorem V.66, and the elementary inequality 
|a + 6|'' < 2'^~^{\a\^ + \b\'') for all a,6 G M that there exists a constant Kk G M+ such that 

i;(r,||^)<i^Ji?(|lyo||^) + PlirilC|l2^t+(2'^l|i?l|^- + Pilf) tEmf,)ds 



Using the Gronwall lemma and arguing analogously to the case k ~ 1 now completes the 
proof. □ 

Proof of Theorem 4.8. Proposition 4.7 ensures the finiteness and local boundedness 
of the first absolute moment needed in the following. From the defining stochastic differ- 
ential equation (3.4), we have 



Yt=Yo+ (Bn_+n_B*)ds+ / A{Ys^+C)^^^d[L,L]l{Ys-+Cy^^A*. 
Jo Jo 

Therefore, 

E{Yt) = E{Yo)+ f {BE{Ys) + E{Ys)B*)ds + aL f AE{Ys + C)A* ds, (6.5) 
Jo Jo 

using a Fubini argument, stochastic continuity, the variant of the compensation formula 
given in Lemma 6.9 and the observation that E{[L,L]\) = J^^xx*i'{dx) = var(Lj) is 
implied by equation (2.1). Thus, 



vec (^E (^J^ A{Ys- + C)i/2 d[L, L]l (n_ + C)^/^A 

= e(^J\a ® A)((y,_ + C)i/2 ® + C)i/2) dvec([L, L]l) 
t 

{A ® A)E{{Ys. + Cfl^ ® (y,_ + C)i/2) vec(i;([L, L\\)) ds 
-GL I {A®A)E{{Y,_+Gf'^®{Y,.+Cf/^^cc{Id))ds 
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^ AE{Y,^ + C)A* ds^ . 



: (Ti vec 



Equation (6.5) therefore implies the following differential equation after vectorizing: 
^E{veciYt)) = SSElYcclYt)) + aL{A® A) vee(C). 

Solving this ODE establishes (i). 

Turning to (ii), the assumed second-order stationarity and (6.5) imply that 

BE{Yq) + E{Yq)B* + aLA{E{YQ) + C)A* = 0. 

The rest is just a matter of rewriting this linear equation. □ 

Proof of Lemma 4.10. From Assumption 4.2, we have that 



(tlII vee(/d)||i3,s 



< I \\YCc{xx*)\\B.,sMdx). (6.6) 

B,S 



For A; = 1, condition (4.4) becomes \\S\\1\\S-^\\IK2.b\\A<»A\\b,s /r^ |1 vcc{xx*)\\b.sM(^x) < 
-2A. Using (6.6), |15||2||5-i||2 > 1 and the fact that i^s.sH vec(/d)|lB,s > ||/d||2 = 1 due 
to the definition of K2,Bi one obtains 

ctl\\A® A\\b,s = adiS-^ ® S-^)iA® A){S ® S)\\2 < -2X. 

Let now be any eigenvalue of ^ and note that {S^^ ® S~^){B ® Id + Id^!^ B){S ^ S) 
is diagonal. Thus, the Bauer-Fike theorem (see [23], Theorem 6.3.2 and its proof, for 
instance) gives that there exists a /i S (j(B (E) Id + Id 'S^ B) such that 

|3fi(^)-5R(/i)| < l/i-Al < \\{S-^<ES-^){^~B®Id-Id®B){S®S)\\ 
= aLWiS-^ <E S-^){A(S) A){S <S) S)\\2 < -2X. 

Hence, 3?(/i) < max{5R(/i) : p. G a(B ® Id + Id^ B)} — 2X = because the maximum equals 
2A due to a{B (g) Id + Id <^ B) = (t{B) + (t{B) and the definition of A = max(5R(CT(S))). 
Therefore, a{^) C (-oo, 0) + iM and d§ is invertible. □ 

Proof of Theorem 4.11. We only prove (i), because the proof of (ii) proceeds along 
the same lines, noting that the finiteness is ensured by Proposition 4.7. 

The equality acovy(-) = acovy(-) is obvious. Due to the second-order stationarity, we 
have 

acovY(/i) = cov ^vec ^Yo + j {BYs- + Y^-B*) ds 

+ £ A{Y,. + C)i/2 d[L, L]l (y,_ + Cy/^A*'^ , vec(yo)) 
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= var(vcc(yo)) + ^ [B^h + h^B) vcc(n_) vec(yo)* dsj 
-E(^j^{B® la + Id®B) vec(y,_) ds^ i;(vec(yo)*) 
+ {A®A){{Ys-+Cf'^®{Ys-+Cf'^)AYCc{[L,L]l) vec(yo)*) 

-i?^\A®A)((y,_ + c)i/2®(i;_ + c)i/2)dvec([i,i],^)^ii;(vcc(ro)*) 

= var(vec(yo))+ / {B ® 1^ + h® B)E{YCc{Y,)YCc{YoY)ds 
Jo 

t 

{B®Id + Id® B)E{vec{Ys)) ds E{veciYo)*) 
OL I (A®A)i;(((n-+C)i/2®(n_+C)i/2)vcc(/d)vcc(Fo)*)ds 

{A(E)A)E{{{Ys^+Cy/^(E){Ys-+Cy/^)vec{Id))dsE{vec{Ya)*) 



(TL / {A® A)E{Yec{Ys+C)vec{Yo)*)As 
Jo 

-<JL I {A®A)E{ycc{Ys+C))E{YCc{Yo)*)ds 



+ var(vec(lo)) + [ (B ® + Id ® B) acovy (s) ds 

^0 

= var(vec(yo)) + / [B ® Id + Id® B + cjlA® A)acoYY{s)ds, 
Jo 

where we have used a Fubini argument, Lemma 6.9 and E{[L,L]i) = (JLid- Regarding 
the use of Lemma 6.9, we observe that \\{Ys-+Cy/^ ® {Y,.+Cy/'% = HXs- +C||2 and 
hence the required local boundcdness is ensured by the second-order stationarity of Y . 

The ordinary differential equation (4.10) is now immediate and to conclude the proof, 
it suffices to note that acovy(O) = var(vec(Fo)) and thus solving the ODE gives 

acovy(/i)=acovy(/i)=e(^®-^'^+^''®^+'^^^®'^)''var(vec(yo)), h>0. □ 

Proof of Lemma 4.14. Let e be as in the definition of type G and let ly^ be its Levy 
measure. Then, by [2], Proposition 3.1, L has Levy density u{x) = /jj^ (j>d{x;Tld)i'e{dT) , 
where </>£;(•; E) denotes the density of the d-dimensional normal distribution with variance 
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E. Hence, 

E{[vec{[L,L]^),vec{[L,L]^)]l) = [ {xx*) ® {xx*)u{x) dx 

= I I (xx*) ® {xx*)(l)dix]Tld) dxiyeidr) 

= / T^i^,{dT){Id2 +Kd + vcc{Id)vcc{Id)*), 

using [30], Theorem 4.3. Now, set := T^i^e(dT) and note that the finiteness follows 
from the definition of type G and the assumed finiteness of the fourth moment of L. □ 

Proof of Theorem 4.15. Proposition 4.7 again ensures the existence and local bound- 
edness of the second moment needed in (i). 

The definition of quadratic variation (see [5], Lemma 5.11, for a special version in the 
context of matrix and vector multiplication) implies that 

vec(rt) vec(yt)* = vec(yo) vec(ro)* 

+ / ((S®/d+/d®S)vec(r,_)vec(y,_)* 
Jo 

+ vcc{Ys-)vcc{Ys-)*{B* (E)Id + Id<^B*)) ds 

+ [\a®A){(C + ®{C + Ys-f'^)dYCc{[L,L\l) vec(K,_)* 

Jo 

+ / Ycc{Y,^)dYcc{[L,L\ly{{C + Ys-f'^®(C + Ys-f^^){A*®A*) 
Jo 

+ [vec(y),vec(r)]t. 
Moreover, setting Vt = ((C + Yt)'^/'^ ® {C + Ft)^/^) we obtain 

[vec(y),vec(y)*]t = / (A ® A)V,_ d([vec([i, vec([L, L]')]°)Vs-(A* ® A*) 
Jo 

= 11 {A® A)Vs^YCc(xx*)\cc(xx*)*Vs-(A* ® A*)pLL(ds,dx), 
Jo Js.'' 

since vec(Yt) is the sum of an absolutely continuous component and a pure jump process 
of finite variation. 

Using a Fubini argument, the stochastic continuity. Lemma 6.9 and the Assumptions 
4.2, 4.3 made on the moments of j^l, we obtain 

i;(vec(yt)vec(yt)*) 
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E{yec{Yo)vec{Yor) 

+ f iiB (g, Id + Id(^ B)E{Yec{Ys)Yec{Y,)*) 
Jo 

+ E{vcc(Ys)vcc{Ysy){B* (E)Id + Id<E) B*)) d.s 
+ aL I [A® A)E{ycc{C + Ys)YCc{Ysy)ds 



+ (TL / £;(vec(n)vec(C + y,)*)(A*® A*)ds 

+ / {A® A)E{Vs-PL{Id- + Kd + Ycc{Id)vcc{Idy)Vs-){A* ® A*)ds. 
Jo 

With the definition of , it follows that 

E{Vs-Id^Vs-) = E{V^) = E{{C + Y,) g,{C + Ys-)) 
= QE{vcc{C + y,) vec(C + n)*), 
E{Vs- vec(/d) vec(/d)*Vs-) = -B(vee(C + Ys) vec{C + F,)*), 
E{Vs-KdVs-) = KdE{{C + Ys) (g{C + Y^)) 

= KdQEivcc{C + Ys) vec(C + Ys)*), 

using [30], Theorem 3.1 (xii), in the last identity. Inserting these formulae into the above 
result and noting that in the stationary case, the integrands need to sum to zero gives 
(4.17). Vectorizing then immediately establishes (4.18). 

Likewise, we obtain (4.16) in the non-stationary case by inserting the formulae above, 
vectorizing and differentiating. □ 

Proof of Lemma 4.17. We have 

PL 1 1 vec ( /<j2 + /\ d + vec (/d ) vec (/d ) * ) 1 1 ^ 

vcc{{xx*) (g) {xx*))iyL{dx) (6.7) 

' B,S 

vec{{xx*) (g {xx*))\\g~'„i^L{dx) = / \\vec{xx*)\\% gi^Lidx) 



< 

since the definition of || • \\]^ implies that 

||vec((a;a;*) g) {xx*))\\^^ = \\S~'^{xg)x) g) {xg)x)\\2 



US-^ g) S-^){x (g x)\\l = li vec(a;x*)|| 
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using the fact that ||z(g) z||2 = \\z\\l for all z G ]R'^^ 
For k = 2, the condition (4.4) becomes 

2\\S\\1\\S-'\\IK2,bU®A\\b,s f \\yec{xx*)\\B^s'^L{dx) 

JR'' 

+ \\S\\t\\S-Y2KlB\\A®A\\%s I ||vec(.Ta;*)||^,5i.i(dz)<-4A. 



Using (6.7) and results from the proof of Lemma 4.10 gives 
2ctl||A® A||B,s + A"2%P®^lis,5Pil|vec(/d2 +A'd + vec(/d)vec(/d)*)|l5;^ <-4A. 

Combining \\{A®A)®Id2 + I^p ®{A® A)\\^g < 2\\A(^A\\b,s and \\A®A®A(»A\\^g = 
\\A(E} AW^ g, which are elementary to prove, with (4.19) leads to — 3§ ® 1^2 ~ 1^2 (g) 

mB:s<-^^- 

Since S~^{^ ® 7^2 + Ij^i ® SS)S is diagonal and max(3fJ(cr(^ ® 1^2 + 1^2 ® SS))) = 4A, 
the Bauer-Fike theorem (see [23], Theorem 6.3.2 and its proof) and arguments as in the 
proof of Lemma 4.10 complete this proof. □ 

Proof of Lemma 4.18. That 5* is unitary implies that K2.B =■ 1 and all the norms used 
are actually the Euclidean norm or the operator norm induced by it. Hence, we have to 
show that 

1 1 Q + /Crf Q + 44 1 1 2 < 1 1 vec(/rf2 + Ki + vec(/d) vec(/d) * ) || 2 • 

For d—\, one calculates both sides to be equal to 3. 

In general, we know from the fact that K,d and Q are permutation matrices that the 
operator norms are 1. Hence, || Q + /C^Q + /d*" II2 < 3. Furthermore, the entries of 
and vcc{Id)ycc{Id)* are either 1 or 0. Therefore, || vec(/d2 + Kd + vec(/(i) vec(/d)*)||2 > 
II vec(/d2)||2 = d. This shows the inequality for d > 3. 

In the remaining case d= 2, we have 

||vcc(/d2 +/<rf + vec(/d)vec(/d)*)||2> ||vec(/d2 + vec(/d) vcc(/rf)*)||2 = \/T2 > 3, 
which again establishes the claimed inequality. □ 

To prove the asymptotic second-order stationarity, we need the following general lemma 
on differential equations which is elementary, but not to be found in the literature, to 
the best of our knowledge. 

Lemma 6.10. Let/:M+->M'' be continuous and A e Md{R) with a (A) C {-00, 0) +iR. 

If limt^oo f{t) = ^ with ^ G M'', then, for any initial value xq G W^, the solution x to 
the differential equation 

dx{t) 



dt 

satisfies limt^oo x{t) — —A^^S^. 



■Axit)+f{t) 
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Proof. It holds that 



Smce 



Hm / e^(*-'*)^ds = -A-i^ and hm e^*xo = 0, 



it suffices to show that hmt_j.oo |1 /q e'^'-* '^Hfi^) ^ O'^^W = 0- Fix e > 0. There exist 
t*,t** > with t* < t** such that \\f{t) - C\\ < s for ah t > t* and 



< e 



for all t > t* 



Hence, 

t 







< 1 



^'^(^-^^llds £< 1 



ds 



Since the last integral is finite and e was arbitrary, this completes the proof. 



□ 



Proof of Theorem 4.20. (i) follows from Proposition 4.7 and Theorems 4.8, 4.11 and 
4.15. 

Regarding (ii), Proposition 4.7 ensures that i?(||yt|p) <oo for all t€ M+. The conver- 
gence of the expectation has already been noted in Remark 4.9 and the convergence of the 
variance follows from (4.16) and the previous lemma. (4.12) then implies the convergence 
of the autocovariance. □ 



6.3. Proofs for Section 5 

Proof of Proposition 5.2. Proposition 4.7 ensures that i?(||Ytj|) and hence £'(||Vt||) 

is finite and locally bounded. Since E{\\Vt''^\\l) = EiWVth), the standard i^.stochastic 
integration theory (see [1], Section 4.2.1, for example) establishes that £'(||G'(|1^) < oo for 
alHe M+. 

If we now let Assumptions 5.1 and 5.2 be satisfied, then (5.2) is clear and (5.5) is a 
straightforward consequence of the ltd isometry. The latter also implies that 

var(Gi)==£;(GiGt) = £;^ V^L^ E{LiL\)Vll^ = {aL + (Tw)E(^j^ Vs-±^ 

= {(TL+(Tw)^E{Vo). □ 

Proof of Proposition 5.4. (i) We first show the finiteness of the moments. Using 
the Euclidean norm and its operator norm, it is clear that i?(||Gt||^) < oo if and only 
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if i^dlGtGj IP) < oo. Denoting the d components of G by Gi with i = 1, . . . ,d, we 
have that £'(||G'f||^) < cxj is equivalent to E{\Gi^t\'^) < oo for all i e {!,..., d}. But the 
Burkholdcr-Davis-Gundy inequalities (see [36], page 222) give that E{\Gi^t\'^) < oo pro- 
vided E{[Gi, Gi]j ) < oo. The latter is, in turn, ensured by [G, G]t|p) < oo simultane- 
ously for all i g { 1 , . . . , d} . 

Next, we observe that [L,L]t = r^t + J^, .Ta;*/iL(ds, dx) with e Sj" and 
that, moreover, [G, G]t = ^^V^i^ A[L,L\sV^i^ or, equivalently, vec([G, G]t) = 
/q V^l^ ®Vl!_^ d\ec{[L,L\s). Since Proposition 4.7 ensures the finiteness and local bound- 
edness oi E{\\Vs\\l) = E{\\Vs (^Vsh) and E{\\Li\\^) the finiteness of [L, ), the 
standard i^-stochastic integration theory immediately gives [G, G]t|p) < oo. 

(ii) If we now let Assumptions 5.1 and 5.2 be satisfied, then (5.6) has already been 
shown in the last proposition. 

It thus remains to establish (5.7). Let e N. The definition of the quadratic 
(co)variation (see [5], Lemma 5.11, in particular) implies that 



JhA 



1/2 



hA 



(h+l}A 
hA 



ih+l)A 
hA 



(h+l}A 
hA 



1/2 



hA 



dL„ 1 dL: V, 



1/2 



with 



{h+l)A 
hA 



(h + l}A 
hA 



vli'dL, 



{h+l)A 



Vli^d[L,L]sV, 



1/2 



hA 



{h+l)A 



V.^-ds 



hA 



(h+l)A 



vli^d[L,L]lV, 



1/2 



We now condition upon Jvv, and denote by Y{y, {Lr — Lt„)r>to,tQ,t) with to e M+ the 
solution of 

dYt = {BYt- + Yt-B*) dt + A{Yt- + Cf'^ d[L, L\l {Yt- + Cf'^A* 

for t>tQ with Yfg = y. Furthermore, we denote by £'L.to(') the expectation taken with 
respect to {Lr — Lta)r>to only. Using Theorem 4.4, E{Li) = and the fact that the 
increments of (Lr)r>A arc independent of J^a, one obtains 



E 



{h+l)A 
hA 



Vli'dL, 



hA 



1/2 



^E 



{h+l)A 
hA 



YlL'dU 



^A 

1/2 



hA 



Ya 



i?L,A( / / V^[l_dLsdL:V^il_ 
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r(/i+l)A j-s 
IhA JhA 

where VY^,t := ^(^^A, (^r — iA)r>A, A, t) + C and, likewise, 

f-(/!.+l)A 



/ /• /!.+! A / /•« \ 

e( ( v^L'dLMLit!-' 

\JhA \JhA / 



-Fa = 0. 



Moreover, using the moment assumptions and the compensation formula, we have 

n(/i+l)A /.(/!.+l)A 

^A 



E[aw I Vs-ds+ I VIL^ d[L, L\ylL^ 

\ JhA JhA 

r{h+l)A 

= {ol+ow) j £'L,A(VVA,s)ds. 
JhA 

Equation (4.7) implies that 

/.{/i+l)A 

/ £;L,A(vec(WA.,))ds 

JhA 

(h+l)A 

(vec(C) + e^('*~'^)(vec(rA) + aL^-^{A A) vec(C))) ds 

HA 

(h+l)A 

aL^-^{A®A)vcc{C) As 

hA 

= A vec(C) + ^-ie^^''(/d2 - e-®'^)(vec(rA) - i^(vec(ro))) + AS(vec(yo)). 
Combining the above results, we get 

i?(vec(G,+iG;;+i)(vec(GiGD)*) 

= £;(ii;(vec(G„+iG,:+i)|J-A)(vec(GiGt))*) 

= £;(vec(GiGD)(i5(vec(GiGD))* + (^l + aw)e-®^'^^-^(/rf2 - e"^^) 
X (£;(vec(rA)(vec(GiGt))*) -i;(vec(ro))(£;(vec(GiGt)))*). 

Using the stationarity of F, this establishes (5.7). □ 
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